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Abstract. The presence of obscuring material (or a dusty 'torus') in active galactic nuclei (AGN) is central to the unification 
model for AGN. Two models, the multi-population model for radio sources and the receding torus model, are capable of 
describing observational properties of powerful radio galaxies and radio quasars. Here, I study the changes of the opening 
angle of the obscuring torus with radio luminosity at 151 MHz, [O III] emission-line luminosity and cosmic epoch aiming to 
discriminate between two working models. An analytical expression relating the half opening angle of the torus to the mean 
projected linear sizes of double radio galaxies and quasars is derived. The sizes of powerful double radio sources taken from the 
combined sample of 3CRR, 6CE and 7CR complete samples are used to estimate the torus opening angle. I found a statistically 
significant correlation between the half opening angle of the torus and [O III] emission-line luminosity. The opening angle 
increases from 20° to 60° with increasing [O III] emission-line luminosity. This correlation is interpreted as direct evidence of 
the receding torus around central engines of powerful double radio sources. 

Key words. Galaxies: active - Radio continuum: galaxies - Methods: analytical 



1. Introduction 

Spectropolarimetric observations of Seyfert 2 galaxies 
(Antonucci & Miller 1985; Miller & Goodrich 1990) sug- 
gested an orientation-based unification scheme for active 
galactic nuclei (AGN). An optically thick molecular torus 
obscuring the AGN permits a unification of Seyfert 1 and 
Seyfert 2 galaxies by proposing that they are intrinsically 
the same objects viewed at different angles to the axis of the 
torus. This model was invoked by Scheuer (1987) and Barthel 
(1989) to account for the orientation-dependent appearance of 
powerful radio sources, which can appear as radio galaxies 
or quasars. According to the standard picture, the torus hides 
the AGN and broad-line region of a radio galaxy, whereas 
for quasars these are viewed directly for a quasar down the 
ionization cone. The opening angle of the torus defines the 
cone angle and the half opening angle is the critical angle 
(6 C ~ 45°; Barthel 1989) at which the transition between 
radio galaxies and quasars occurs. The narrow-line regions 
have extended kpc-scale structure and therefore are visible for 
any orientation of the torus. There is considerable multi-band 
observational evidence supporting the unification scheme 
for radio galaxies and quasars. This comes from the X-ray 
detection of inverse Compton scattered photons (Brunetti et 
al. 1997) and polarized (scattered) emission lines from hidden 
quasars in some radio galaxies (Antonucci & Miller 1985), 



from the direct detection of unabsorbed soft and hard X-ray 
emissions from quasars (e.g. 3C 109, Allen and Fabian 1992; 
Cyg A, Ueno et al. 1994; 3C 295, Brunetti et al. 2001; 3C 
265, Bondi et al. 2004), from the detection of comparable 
far-infrared emission from the high-redshift (z > 0.8) radio 
galaxies and quasars (e.g. Meisenheimer et al. 2001; Haas 
et al. 2004), from the difference of the mid-infrared spectral 
energy distribution in the broad-line and narrow-line AGN 
(Siebenmorgen et al. 2004), from the detection of superluminal 
radio knots on parsec-scales (Vermuelen & Cohen 1994), and 
from the presence of supermassive black holes and relativistic 
jets in both radio galaxies and quasars (McLure & Dunlop 
2002; Laing 1988; Garrington et al. 1988). 

A single opening angle can not account for two pieces of 
observational evidence: (i) at a given radio power the [O III] 
lines are weaker in radio galaxies than in quasars (Lawrence 
1991; Simpson 2003), and (ii) the increase of the quasar frac- 
tion with increasing radio luminosity (Hill et al. 1996; Willott 
et al. 2000; Grimes et al. 2004). Lawrence (1991) suggested 
the receding torus model in which the dust evaporation radius 
of the torus depends on the luminosity of the photoionizing ra- 
diation of the AGN, which is assumed to be independent of the 
height of the torus. In this scenario the opening angle of the 
torus increases with increasing luminosity of the ionizing ra- 
diation, Lphot, from the AGN. This model accounts for the fact 
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that quasars are brighter in L p h t than radio galaxies (Lawrence 
1991; Simpson 2003) in samples selected at low radio frequen- 
cies, and that the fraction of quasars increases with radio lumi- 
nosity. Gopal-Krishna et al. (1996) showed that a wider open- 
ing angle with increasing radio luminosity could explain the 
problem of unification posed by the difference between the ra- 
dio luminosity-size correlations for radio galaxies and quasars. 
In addition, Willott et al. (2000) found that the fraction of 
quasars decreases at low luminosities. They argued that this 
can be explained either by the receding torus model (decrease 
of 9 C with decreasing radio luminosity), or by the existence of 
a second population of low-luminosity radio galaxies. Grimes 
et al. (2004) used a combined complete sample at 151 MHz to 
show that, in the luminosity-dependent unified scheme, a two- 
population model may explain the increase of the quasar frac- 
tion with Lphot as well as the emission-line differences between 
radio galaxies and quasars. They conclude that the effect of the 
receding torus may be important but it is not clear yet whether 
the receding torus is present in both populations. 

The opening angle of the torus, 9 C , defines the geometry 
of the torus (the inner radius and height) and appears to be 
an important parameter for understanding the evolutionary be- 
haviour of the obscuring material and testing the torus model. 
In the Appendix, I derive an equation that allows the mean half 
opening angle of the torus 9 C to be estimated by means of the 
projected linear sizes of radio galaxies and quasars. In section 
2 the equations for estimating 9 C are presented. The combined 
complete sample of radio sources is described in section 3. In 
section 4, the relations between the critical angle with the 151 
MHz radio luminosity L151, [O III] 5007-A emission-line lumi- 
nosity L[oni] and cosmic epoch are investigated. In section 5, 
the receding torus model is discussed. 

Throughout the paper a flat cosmology (£? m + — 1) with 
non-zero lambda, Q m = 0.3 and Hq = 70 km s Mpc~' is 
used. 

2. Equations for the opening angle of the torus 

Powerful FRII (Fanaroff & Riley 1974) radio galaxies and 
quasars are associated with bipolar relativistic radio jets emerg- 
ing from the AGN which form extended radio lobes on both 
sides of the AGN. The radio jets are almost aligned with the 
axes of the tori in FRII (Cyg A; Canalizo et al. 2004) and FRI 
radio sources (3C270, Jaffe et al. 1993; Verdoes Kleijn et al. 
2001). The half opening angle of the torus is defined by the an- 
gle between the radio axis and the direction at which the divi- 
sion between radio galaxies and quasars occurs. Then the half 
opening angle of the torus can be estimated by means of the 
projected linear sizes of FRII radio galaxies and quasars; lin- 
ear sizes of quasars appear to be systematically smaller than 
the linear sizes of radio galaxies (Barthel 1989) as a result of 
projection effects. In the Appendix, I derive an equation (7) 
for estimating the mean critical angle # cs from the ratio of the 
mean projected linear sizes of quasars Iq and the whole sample 
I. The problem is solved for the sample of FRII radio sources 
having an isotropic distribution of radio axes over the sky. Low- 
frequency radio samples of radio sources are thought to be free 
from orientation biases and can be used in this analysis. 



Table 1. The Kolmogorov-Smirnov test is used to estimate the 
significance level of the hypothesis that the [OIII] emission- 
line luminosity (and radio luminosity) distributions of HEGs 
and quasars (QSOs), HEGs and WQs, and, QSOs and WQs are 
drawn from the same parent distribution. 





K-S test 


^151 


K-S test 


distribution 


(%) 


distribution 


(%) 


QSO-HEG 


99.99 


QSO-HEG 


99.66 


QSO-WQ 


98.77 


QSO-WQ 


99.81 


HEG-WQ 


16.60 


HEG-WQ 


79.47 



Another independent way to estimate the half opening an- 
gle is to consider the number of quasars Nq and the number 
of radio galaxies and quasars Nq+q in the low-frequency radio 
samples, 

rf = arccos(l--i). (1) 

If the unified schemes for FRII radio galaxies and quasars is 
valid then one should expect a correlation between 9 C { and 9 CS , 
which are independent estimates. 

3. The sample and selection criteria 

I use a combined sample of three complete low-frequency 
samples (Grimes et al. 2004) made up of the 3CRR, 6CE 
and 7CRS samples (Laing, Riley & Longair 1983; Rawlings, 
Eales & Lacy 2001; Willott et al. 2002). This contains 302 ra- 
dio sources all having extended radio structure on kiloparsec- 
scales. I then excluded 3C48, 3C287 and 3C343 because they 
are complex with no clear double lobe structure (Laing et al. 
1983). Three other sources were excluded: 3C231 because its 
radio emission is due to a starburst, 3C345 and 3C454.3 be- 
cause their inclusion in the 3CRR sample is a result of Doppler 
boosting of the jet. 

I adopted the classification of radio sources used by Grimes 
et al. (2004): high- and low-excitation narrow-line radio galax- 
ies (HEGs and LEGs), weak quasars (WQs) and broad-line 
radio quasars. The category of weak quasars includes HEGs 
which are objects with weakly or heavily obscured broad-line 
nuclei seen indirectly (e.g., in broad wings of Ha and Paschen 
lines, in optically-polarized light) and unobscured objects with 
weak broad-line optical nuclei. 

There are 39 low-excitation radio galaxies in the sam- 
ple which are believed to be a separate population indepen- 
dent from the high-excitation radio galaxies and radio quasars 
(Laing et al. 1994). Exclusion of these sources and FRI radio 
sources leaves a sample consisting of 237 FRII radio sources. 
The projected linear sizes (the apparent distance between radio 
lobes) of 3CRR sources are adopted from the 3CRR database 2 , 

1 The electronic version of the combined sample is available at 
http://www-astro.physics.ox.ac.uk/~sr/grimes.html 

The electronic version of the 3CRR database can be found at 
http://www.3crr.dyndns.org/cgi/database 
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Fig. 1. The distribution of [OIII] emission-line luminosity 
for 134 high-excitation narrow-line galaxies (open area), 52 
quasars (grey area) and 20 weak quasars (hatched area) in the 
combined 3CRR, 6CE and 7CRS samples. 
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Fig. 2. The distribution of radio luminosity at 1 5 1 MHz for 1 34 
HEGs (open area), 52 quasars (grey area) and 20 WQs (hatched 
area) in the combined sample. 

and the sizes of 6CE and 7CRS sources are taken from the elec- 
tronic version of the combined database 1 . The projected linear 
sizes vary over (1 to 2000) kpc over the range (1.6 x 10 24 to 
4.6 x 10 2S ) W Hz 1 sr~' of the 151 MHz radio luminosity. 

The structure of FRII radio sources is influenced by en- 
vironmental asymmetry, which becomes more important on 
small scales (Arshakian & Longair 2000). I excluded the com- 
pact steep-spectrum radio sources with typical sizes < 20 kpc to 
avoid contamination of projection effects by interaction effects 
between a radio jet and interstellar medium (Barthel 1989). 
Their inclusion or exclusion does not introduce perceptable 
change in the results described below. The final sample con- 
sists of 206 FRII radio sources: 134 HEGs, 20 weak quasars 
and 52 quasars. 

4. Relations for the torus opening angle 

The distribution of [OIII] emission-line luminosities and ra- 
dio luminosities of 206 FRII radio sources is shown in Figs.[0 



35 35.5 36 36.5 37 
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Fig. 3. Mean projected linear sizes of 206 FRII radio sources 
are calculated for three equal subsamples binned by [O III] 
emission-line luminosity. Quasars, weak quasars and HEGs are 
denoted by 'Q', 'W and 'H' respectively. 



and |3 There are 10 radio galaxies and two WQs with power 
< 2 x 10 25 W Hz 1 sr , typical for the transition region for 
FRI and FRII sources. Quasars appear to be more luminous 
than HEGs and WQs in both [OIII] emission-lines and radio. 
The Kolmogorov-Smirnov test (Table ^ shows that the lumi- 
nosity distributions (L[om] and L151) of quasars are different 
from the distributions of HEGs and WQs at high confidence 
level (mainly because of the lack of weak quasars in the sam- 
ple), whilst the L[omj and Z451 distributions of HEGs and WQs 
appear to be drawn from the same population. The later result is 
supported by the linear-size statistics of FRII sources: the mean 
projected linear size of weak quasars, (330 + 50) kpc, is compa- 
rable with the mean linear size of narrow-line HEGs, (320 + 30) 
kpc, and is twice as large as the mean size of quasars (162 ± 19) 
kpc. This relation holds over the entire [O III] emission-line 
luminosity range (Fig. The values of mean linear sizes of 
weak quasars are larger than those of quasars in three equal 
subsamples divided by [O III] emission-line luminosity. The 
same result is obtained when the sample is binned by 15 1 MHz 
radio luminosity. This supports the idea that HEGs and weak 
quasars are the same objects, where the later ones appear to 
have broad emission lines as a result of scattered light from 
the nucleus. Therefore in this analysis, the weak quasars are 
grouped with the HEGs rather than quasars. On the other hand, 
some of the WQs appear to be transitional objects (between 
quasars and HEGs) which are oriented near to the critical an- 
gle where the broad-line region is reddened but not totally ob- 
scured (Laing et al. 1994; Simpson et al. 1999) and hence they 
can be grouped with quasars. Therefore it is important to see 
how an exclusion of 20 WQs influences the final results and 
this will be discussed later in this section. 

For the whole sample, I calculate the number of all radio 
sources Ng+q = 206 and quasars Nq = 52, and their mean 
linear sizes I = (279 ±21) kpc and Iq = (162 + 19) kpc re- 
spectively. Then the critical viewing angles, CS = 42° + 7° and 
9 C f = 42° ± 12°, are estimated using equations (7) and (1). Two 
independent equations give similar results with the error being 
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Table 2. Estimates of critical angles, 6* cs and 8 c f and their errors in equal subsamples binned by L[onrj- The first line of the table 
refers to the two sets of binned data, and the second line shows results of three sets of binned data. 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 

N bm log(Z.[orrrj/W) log(L[onrj/W) N a+Q N Q e cf ±cr Bc[ Fisher test l±a-j 1q±o-- 1q cs ± cr 9a Fisher test 

bin width (°) (%) (kpc) (kpc) (°) (%) 

~~ 2 32.43 -35.81 35.51 103 12 28 ± 35 48^6 341 ± 36 111 ±25 22 ± 6 99A 

35.81 -37.37 37.08 103 40 52 ±12 217 ± 19 177 ± 22 63 ± 14 

"1 32.43 - 35.38 35.08 68 6 24 ± 57 17l 311 ± 11 79 ± 19 17 ±5 9TA 

35.38 -36.13 35.90 68 15 39 ± 24 308 ± 43 175 ± 39 41 ± 12 

36.13 - 37.37 37.09 70 31 56 ± 13 225 ± 23 178 ± 25 59 ± 14 



Columns: (1) - The number of bins; (2) - [O III] luminosity bin width; (3) - Logarithm of the central value of the L[omj in the bin; (4,5) - 
Number of radio galaxies and quasars in the log L[ iirj bin; (6) - The half opening angle of the torus estimated by fraction of quasars (eq. 
1), and its standard error cr Scf is calculated assuming that cr Na ^ Q = -^Ng+q and ct Nq = -^/Nq; (7) - The confidence level of rejecting the null 
hypothesis that half opening angles # c f of the torus estimated in different L[om] bins are equal. The variance analysis (the Fisher test) is 
used with two degrees of freedom, N Vm - 1 and N biu (N G+ Q - 1); (8,9) - The mean linear sizes of all sources (radio galaxies and quasars) 
and quasars, and their standard errors; (10) - The half opening angle of the torus and its standard error from eqs, (7) and (8); (11) - The 
confidence level for rejecting the null hypothesis that half opening angles 9 CS of the torus are equal in different Z40111] bins. 
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Fig. 4. The torus half opening angle calculated in three equal 
subsamples binned by [O III] emission-line luminosity (see 
Table|2j- Critical angles 8 CS estimated by the mean linear sizes 
of radio galaxies and quasars (eq. 7) are marked by filled cir- 
cles, and open circles denote the critical angles 9 c f calculated 
by the fraction of quasars (eq. 1). A Fisher test rejects the null 
hypothesis that CS are equal at 97.4 % confidence level. The 
filled circles are offset horizontally for illustrative purposes. l<x 
error bars are presented. 

smaller for the mean linear size statistics (see Appendix for ex- 
planation). The true value of 8 CS lies between ~ 20° and ~ 60° 
at 3cr (99.7 %) confidence level. There is a spread in the values 
of the critical viewing angle, which indicates that the simple 
unification scheme with constant 8 C is not likely - there is some 
distribution of C . It would be interesting to see if the spread in 
8 C correlates with any other physical property of radio sources. 

To examine the 6 C — L[oni] dependence, I divided the 
sample into two equal subsamples (118 and 119 sources) by 
log(L[om]/W) = 35.83 and calculated 8 CS and G c f for each sub- 
sample (Table 0. Both, 8 CS and 8 c f, increase by a factor of two 



Table 3. The significance (Fisher test) of the null hypothe- 
sis that the 6> cs are equal in two equal subsamples divided by 
£[Oiii]> ^451 an d redshift z, respectively. Probabilities are cal- 
culated for: (1) the sample which includes 20 WQs grouped 
with HEGs rather than quasars, (2) the sample consisting only 
of HEGs and quasars, and (3) the sample which excludes 12 
transitional sources with log(Li5i /W Hz -1 sr~') < 25.3. 



Relation 


Fisher test (%) 


Fisher test (%) 


Fisher test (%) 


plane 


(1) 


(2) 


(3) 


#cs ~ Aoni] 


99.4 


99.3 


98.6 


Ocs ~ ^151 


63.6 


71.7 


60.2 


cs -z 


4.8 


32.5 


38.7 



times from low to high L[onr|. The Fisher test is used to test 
the null hypothesis that the mean critical angles are equal at low 
and high L[onrj. For # cs the null hypothesis is rejected at high 
confidence level (99 %), though it can not be rejected for 6> c f (49 
%) mainly because of the large errors involved (see Appendix). 
I also studied the 8 C - Z40111] plane for three equal subsamples 
(Tabled Fig.0}. The mean half opening angle increases from 
~ 20° to ~ 40° and to ~ 60° in three subsamples binned by [O 
III] emission-line luminosity, and the difference between the 
angles is still very significant (97 %). The difference becomes 
less significant (< 85 %) when the sample is divided on four or 
more subsamples. The division into subsamples smaller than 
~ 50 objects does not allow a statistically significant correla- 
tion between 6 C and L[oni] to be investigated. Also, the small 
samples may introduce fluctuations of 9 C due to statistical fluc- 
tuations in the orientation of radio sources (Saikia & Kulkarni 
1994). 

I repeated the analysis for the 8 C - L151 and 8 C - z relations. 
Both CS and 8 c f estimated for two equal subsamples binned by 
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Fig. 5. The half opening angle of the torus calculated in three 
equal subsamples binned by radio luminosity. Designations are 
the same as in Fig. [4] The null hypothesis of equal critical an- 
gles CS (filled circles) can not be rejected (66%). 



Li5i increase from low- to high-luminosities. This result was 
expected because the relation L[onrj i[ 5 j 4 holds for the com- 
bined 3CRR, 6CE and 7CRS samples (Grimes et al. 2004). A 
Fisher test shows that the increase of 8 CS with Ligi is not sta- 
tistically significant (Tabled Fig- 13- No significant evolution 
of 9 CS or 6> c f is found with cosmic epoch z. 

The results of this section stand even after excluding the 
WQs from the sample (Tabled). This indicates that the positive 
correlation in the 6 CS - L[oni] plane is mainly due to HEGs and 
quasars. The significance of the correlation is still high after 
exclusion of transitional sources from the sample of 206 radio 
sources. 

Let us consider the selection effects which might cause the 
positive correlation 9 C - L[oni] seen in Fig.|4] Although there 
are problems with identifying the broad-line radio galaxies as 
low-luminosity quasars, Hardcastle et al. (1998) suggest that 
some of the low-luminosity broad-line radio galaxies are true 
quasars which are not sufficiently bright to be classified as such 
in the optical. If so, then reassigning them to the HEGs rather 
than to the quasars will decrease the fraction of quasars, re- 
sulting in small opening angles and the consequence of this is 
that 6 c f = 24° is underestimated in the Ljoni] low-luminosity 
bin (Table |3- In a L[oni] high-luminosity bin the fraction of 
quasars may increase if the [O III] emitting region is partially 
obscured in radio galaxies causing a factor of 5 to 10 depression 
of the [O III] emission (Hes et al. 1996) which may result in the 
overestimation of 6 C {. The selection effects at low- and high- 
luminosities may reproduce a positive correlation between 9 c f 



and Li 



[OIII] 



To understand the bias introduced in the fraction 



of quasars one may consider the independent measurements of 
the opening angle 8 CS (Table|2). The fact that 9 CS is independent 
of the fraction of quasars (eq. 7), 9 CS positively correlates with 
L[om] (Fig. and 6 CS =s c f across Lpiii] range indicates that 
the 8 C f - L[oni] correlation is real and if there is a bias in the 
fraction of quasars it is not significant. The x 2 test is used to 
estimate the significance of the null hypothesis that two inde- 
pendent measurements of 9 CS and 9 c f are equal across the L[oni] 



and Li5i ranges. For two and three sets of binned data the val- 
ues of 9 CS and 9 C { are not significantly different (P < 80 %). 
Figs.|4]and |5]show that 9 CS =* cf in different bins of L[oni] and 
Z451. The equality of two independent estimates confirms the 
validity of the orientation-based unification scheme for radio 
galaxies and quasars and it indicates that the relative mean size 
of quasars (Iq/1) and the relative number of quasars (Nq/Ng+q) 
are indeed good measures of the critical viewing angle in the 
low-frequency radio samples. 

Another possible selection effect is related to the measured 
projected sizes of FRII sources which may be overestimated as 
a result of projection of the volume of radio lobes. It is clear 
that independent of the plausible geometry of radio lobes the 
projection effects may be important for sources having small 
inclinations of radio axes to the line of sight. The fact that the 
exclusion of small sizes (< 20 kpc, see section 3) does not 
change the final results indicates that the projection effect is 
not likely to influence strongly the estimates of # cs . 

The assumption that the axis of a dusty torus is aligned with 
the radio axis of FRII sources can be partially relaxed (e.g., 
Capetti & Celotti 1999) in real life allowing some distribution 
of angles between the axes of the torus and the jet. In this sce- 
nario there should be a difference between the true critical an- 
gle (# ct ) defined by the geometry of the torus and the critical 
angle CS (eq. 7) estimated by means of projected sizes of FRII 
sources. It is important to understand how this difference varies 
depending on the value of true critical angle and whether it can 
lead to a positive correlation found between # cs and Lrpm]- To 
test this, I model random angles between the axes of the torus 
and the jet by generating a pair of 5000 random directions: for 
every random axis of the torus the random jet direction is gen- 
erated within the specified cone angle 2<5 max . For simplicity, a 
single linear size (300 kpc) for radio galaxies and quasars is 
adopted. Then, given the true critical angle (9 ct ), I generate the 
projected linear sizes for both quasars and all sources using the 
inclination angles of the jet. The mean sizes of quasars and the 
whole sources are used to estimate the critical angle 9 CS (eq. 7). 
The deviation of estimated critical angle 9 CS from the true 9 ct 
one is presented in Table |4] for different (5 max . <5 max = corre- 
sponds to the case where the axes are aligned and here there is 
an excellent agreement between 9 ct and 9 CS , which provides a 
good test of the formalism involved in equation (7). For larger 
values of 5 max the estimated critical angle appears to be larger 
than the true critical angle (9 CS > 9 ct ) (Table@}. One may con- 
clude that the effect of misalignment of axes of the torus and 
the jet leads to the overestimation of critical angles 9 CS , which 
is higher for small ct . If there is a significant misalignment of 
axes then the critical angles in the 9 CS - Lpiii] relation plane are 
overestimated with the effect of making flatter the slope of the 
real 9 ct - Z40111] relation. This demonstrates that the misalign- 
ment effect makes the real correlation weaker and thus can not 
produce the positive correlation found between 9 CS and L[oni]- 

5. Discussion and conclusions 

The key result of the previous section is that the opening an- 
gle of the ionizing cone of the torus becomes larger at high [O 
III] emission-line luminosities. The half opening angle tan 9 C = 
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Table 4. The values of true critical angles 8 ct and their deviations (0 CS ) simulated for different half cone angles <5 max within which 
the angles between the axes of the torus and the jet are distributed. 
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= 0° 


Omax 
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SmiK = 20° 


<^max 


30° 
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20 


19.8 


±0.4 


21.6 ± 


0.5 


26.5 ±0.7 


34.9 ± 


0.9 


40 


39.5 


±0.5 


40.2 ± 


0.5 


42.8 ±0.6 


47.0 ± 


0.7 


60 


59.3 


±0.6 


59.9 ± 


0.3 


61.2 ±0.7 


63.4 ± 


0.7 



2r/h depends on the inner radius, r, of the torus and its height, 
h (the thickness of the dusty torus). In extreme cases, the posi- 
tive correlation between 9 C and Z^oiirj can be interpreted as (i) 
a systematic increase of the radius of the torus with L[orrn, or 
(ii) decrease of the height with L[oni] ■ The first relation natu- 
rally follows from the receding torus model (Lawrence 1991; 
Hill et al. 1996). In this model (i) the ionizing radiation from 
the AGN evaporates the circumnuclear dust forming the inner 
wall of the torus at a distance r which increases as r ~ Lr^j^ 
(where Z40111] is assumed to be a measure of the ionizing lu- 
minosity), and (ii) the height of the torus h is independent of 
L[0 in] ■ There is observational evidence in favour of an r-L[o m] 
relation. Minezaki et. al. (2004) used multi-epoch observations 
of the Seyfert 1 galaxy NGC 4151 to show that the time lag 
between the optical and near-infrared light curves grows with 
optical luminosity as logAf oc L 5 . They interpreted this as 
thermal dust reverberation in an AGN to relate a At to the inner 
radius of the dusty torus (Af oc r). Willott et al. (2002) found 
that 3C/6C quasars have higher submillimetre luminosities by 
a factor > 2 than radio galaxies. They argued that this is in 
quantitative agreement with the receding torus model if there 
is a positive correlation between optical and far-infrared lumi- 
nosities. Another piece of evidence comes from the studies of 
the mid-infrared spectra of 3C radio sources (Freudling et al. 
2003). They showed that the bands of polycyclic aromatic hy- 
drocarbons (PAHs) are weak in broad-line radio galaxies and 
they are much stronger in narrow-line radio galaxies of simi- 
lar luminosities. They interpreted this difference as a result of 
the heating radiation from the central nucleus: the hard radia- 
tion destroys PAHs in the broad-line regions and heats larger 
dust grains at intermediate distances. The spectra of broad-line 
galaxies and quasars originate in the broad-line region where 
only few PAHs survive. In the narrow-line radio galaxies this 
region is obscured and the spectra are dominated by cooler dust 
at larger distances where more PAHs can survive. Simpson & 
Rawlings (2000) showed that the near-infrared spectral slopes 
of 3CR quasars are steeper for less luminous objects and that 
the fraction of moderately obscured, red quasars decreases with 
increasing radio luminosity. Both relations are in agreement 
with the receding torus model. 

As indicated by Simpson (2003) there are several lines of 
evidence indicating that the height of the torus is not a strong 
function of luminosity. If so, the half opening angle of the torus 
tan# c = 2r/h oc t^/A should increase with [O III] emission 
line luminosity even if there is a spread in h. This trend is in 
agreement with the positive correlation found between the half 



opening angle of the torus and [O III] emission-line luminosity 
(Fig-EJ- To fit the 9 C - L[oni] positive correlation, I used the 
relation predicted by the receding torus model, tan 6 C oc Lf om ^, 

where A is a free parameter. Using the measurements of 6 CS and 
9 c f, and L[om] (see Table |3 I calculate the power-law indices 
A s = 0.35 + 0.08 and A f = 0.26 + 0.02 corresponding to 8 CS 
and 9 c f respectively. The value 0.35 is in agreement with the 
predicted value of 0.5 (r oc i^oun) within ~ 2<x significance 
level. It is likely that the use of the averaged critical angles and 
small number of points result in a smoothing of the real rela- 
tion, which results in smaller values of A s and Af compared 
with the predicted one, 0.5. Larger values of A s = 0.44 and 
Af = 0.27 are found when the number of Z40111] bins are dou- 
bled. The range of A ~(0.3 to 0.4) is a plausible estimate for 
both A s and Af . If these values are correct then in order to fit 
the receding torus model one should expect that the height of 
the torus also correlates positively with the emission-line lu- 
minosity in the form h oc Lp UI] where B ~(0.1 to 0.2). From 
the consideration of the ratio of the near-infrared luminosity 
to the bolometric luminosity for Palomar-Green quasars, Cao 
(2005) found significant correlation between the torus thick- 
ness and the bolometric luminosity, h oc L^ 7 * 05 , which is 
somewhat stronger than is expected from my results. A larger 
sample of radio sources is needed to confirm the increase of the 
torus thickness with photoionizing radiation of the AGN. 

The second relation, i.e. the negative correlation between 
h and L[oni] which may be a cause for an increase of 6 C 
with L[onrj, seems to be unrealistic since it requires the inner 
radius of the torus to be independent of Lfoirrj. One may 
conclude that the growth of the opening angle of the torus with 
increasing Z40111] is a direct evidence of the receding-torus-like 
structures around central nuclei of FRII radio sources. This 
result apparently breaks the degeneracy (multi-population or 
receding torus) in the interpretation of data in favor of the 
receding torus model. 

The main conclusions of this analysis are as follows: 

1 . On the basis of an inverse problem approach, an equation 
for estimating the opening angle of the torus by means of 
the mean projected linear sizes of FRII radio galaxies and 
quasars has been derived. The small errors involved in esti- 
mation of the opening angle makes the linear size statistics 
a powerful tool for investigating the orientation-dependent 
structures in AGN. 
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2. An application of this equation to the combined complete 
sample of FRII radio sources revealed the following: 

- The mean half opening angle of the torus is 42° + 7° for 
the whole sample. 

- Statistically significant positive correlation is found 
between the opening angle of the torus and [O III] 
emission-line luminosity. The mean half opening angle 
grows from ~ 20° to ~ 60° between L [0 m] ~ 10 35 W 
to 10 37 W. This dependence favours the receding torus 
model and is fitted with the tan 9 C oc L^ u] relation. 

- No significant changes of the opening angle are found 
either with radio luminosity or cosmic epoch. 

- Two independent measurements of the half opening an- 
gle are correlated across [O III] emission-line and radio 
luminosity ranges. This supports the validity of the uni- 
fied scheme of AGN and obscuring torus model. It also 
confirms that the mean relative size and relative number 
of quasars in low-frequency radio samples can serve as 
good estimators of the opening angle of an obscuring 
material. 
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APPENDIX: Equation for estimating the critical 
viewing angle 

Suppose, we have a population of double radio sources with ra- 
dio axes oriented randomly in space. The objects are classified 
as radio galaxies if the viewing angle between the axis of the 
obscuring 'torus' and the line of the sight is larger than a criti- 
cal angle defined by a half opening angle of the torus, 9 > 9 C , 
otherwise, if 9 < 9 C , they are viewed as radio quasars. There 
is observational evidence that the axis of a radio jet is aligned 
with the axis of the torus (Jaffe et al. 1993; Verdoes Kleijn et al. 
2001 ; Canalizo et al. 2003). Assume that radio axes of FRII ra- 
dio sources are aligned with the axes of obscuring tori. Then the 
half opening angle of the torus is measured by the angle (9 CS ) 
between the radio axis and the critical angle 9 C , i.e. 6> c w # cs . 
The task is to derive an equation which allows 9 CS to be esti- 
mated given the distribution function of projected linear sizes 
of radio galaxies and quasars, f(l). 

The projected linear size / = losin9, where Iq e [0,l m ] is 
the deprojected linear size (l m being the largest linear size in the 
sample), and 9 e [0, n/2] is the angle between the radio axis (or 
axis of the torus) and the line of sight. Assume that torus axes 
are distributed isotropically over the sky, K(9) = sin 9, then, 

1 fc^-p 
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where F(lo) is the distribution function of linear sizes. This in- 
tegral equation was first derived by Chandrasekhar & Munch 
(1950) for the problem of reconstruction of the distribution of 
stellar rotation velocities. 

The mean linear size of the whole sample is given by, 



- 4- 

h = -l- 

For quasars (8 < CS ) the equation (2) can be written as, 



fd) = C Q l 



r F(k)di 

£ a k Jll - I 2 



(I < Z m sin 6> cs ), 



(3) 



(4) 



where Cq = 1/(1 - cos 9 CS ). The equation for the moments is, 

; m sin9 cs 

r i' ... 

V n = Cq V 0n , (5) 

J m+1 



71+1 Ii2 _ 12 

o a/ o 



and the first moment, which corresponds to the mean linear size 
of quasars, 



(Io)q - 2/q 



1 - cos ( 



"a 



sin t/ cs cos ( 



(6) 



The radio galaxies and quasars have the same distribution of 
linear sizes, hence (Zo)q = Iq- From (3) and (6) one may recover 
the equation for estimating the critical viewing angle, 



x Iq #cs - sin 6 CS cos ( 
2 ] 1 - cos 8 CS 



(7) 



Given the standard errors of the mean projected linear sizes of 
quasars crj Q and the whole population crj, one can estimate the 
standard error of the critical angle, 



'd6a 

\ dl 



88 Ci 



where, 

<9#cs _ 2 (8 CS - sin6> cs cosfl cs ) 2 
dl * AIq 

<9# cs n (1 - cos0 cs ) 2 



A = sin # cs (0 CS - sin 6 C cos 8 CS ) - 2 sin 8 CS (1 - cos 8 C ). 



(8) 

(9) 

(10) 
(11) 



If there is a spread in 8 CS then eqs. (7-11) can be used for 
estimating the mean half opening angle of the torus 8 CS and its 
standard error cr-g . 

Let us investigate the errors associated with critical angles 
(£> cs and 8 c f) defined by the mean size and fraction of quasars 
respectively (eqs. 7 and 1). For this, I adopt the mean linear 
size of A^g+q = 206 radio galaxies and quasars, I = (279 + 21) 
kpc (see section 4), and the mean size of A^q = 52 quasars, 
Iq = (162 ± 19) kpc, from the combined sample of FRII radio 
sources (see section 3). I assume that the ratios crj/l = 0.075 
and o"j q /^q = 0.099 are constant for any critical angle # cs . 




Fig. 6. Functional relation between the half opening angle of 
the torus and its standard error is presented for equation (8) 
(solid line) and equation (12) (dotted line). 

For any given 8 CS , I calculate crg a from eqs. (8-11) using the 
estimates of Iq (eq. 7) and its standard error (provided that 
their inverse ratio is equal to 0.099). The <Tg cs - 6> cs depen- 
dence is shown in Fig. |6] (solid line). The relation between 
<Tg c[ and 8 c f (dashed line) is drawn on the assumptions that (i) 
o-«c + q = V-^g+q - 14 and ct Nq = ^Nq - 7, and that (ii) the 
ratios <Tnq +q /Ng+q = 0.068 and ctn q /Nq = 0.13 are unchanged 
over 8 c f. Given the 8 c f, one may calculate, 
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■ 

where, 

d8 cf 
dN c+Q 

d0 c f _ 1 
8Nq N g +q 



"G+Q 



d8« 
3Nc 



Nq 



N 2 

G+Q 



1-1- 



2\ 



Nq 
Ng + qI j 



-1/2 
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Nq 

N g+ q 



2\ 



-1/2 



(12) 



(13) 



(14) 



One can see (Fig.|6} that the errors, <Tg cf and crg cs , are compa- 
rable at angles > 50°. At smaller angles, the cr flcr grows expo- 
nentially, reaching the value of 50° at 8 C = 20°, whilst cr e . s 
decreases gradually to zero degrees. The important implica- 
tion of this analysis is that the errors associated with 8 CS are 
small indicating that the linear size statistics is a powerful tool 
for studying the correlations between the opening angle of the 
torus and physical characteristics of AGN in the low-frequency 
radio samples. This approach also can be used to estimate the 
opening angle of structures (around the central engines of dou- 
ble radio sources) showing the orientation-dependent physi- 
cal/spectral characteristics. 



